Dynamic Colloidal Stabilization by Nanoparticle Halos 
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We explore the conditions under which colloids can be stabilized by the addition of smaller parti- 
cles. The largest repulsive barriers between colloids occur when the added particles repel each other 
with soft interactions, leading to an accumulation near the colloid surfaces. At lower densities these 
layers of mobile particles (nanoparticle halos) result in stabilization, but when too many are added, 
the interactions become attractive again. We systematically study these effects -accumulation re- 
pulsion, re-entrant attraction, and bridging - by accurate integral equation techniques. 

PACS numbers: 82.70Dd,83.80.Hj 



Colloidal dispersions - solid particles with radii rang- 
ing from a few nm to a few jj,m, suspended in a liquid 
solvent - are common in nature and widely used in indus- 
try. Blood, paint, ink and cement are typical examples. 
Because interatomic dispersion forces induce effective van 
der Waals interactions with large attractive values at con- 
tact, colloids will irreversibly aggregate, which is usually 
undesirable, unless their surfaces are prevented from ap- 
proaching too closely. The two most common ways to 
achieve this are called steric and charge stabilization^. 
Popular strategies for steric stabilization usually involve 
grafting a layer of polymers onto the colloid surface, re- 
sulting in dense repulsive brushes that prevent close con- 
tact. For charge stabilization, the route most common in 
nature, the colloids have surface charges of the same sign, 
leading to a double layer of microscopic co- and counter- 
ions. Adding this effective repulsion to the intrinsic van 
der Waals attraction results in the famous Derjaguin 
Landau Verwey Overbeek (DLVO) potential^], with a 
metastable free-energy barrier preventing aggregation. 

In an important recent development, a third strategy 
for colloidal stabilization, termed nanoparticle haloing, 
was introduced by Lewis and co-workers 0- By adding 
charged hydrous zirconia nanoparticles of average radius 
3nm to a suspension of (marginally charged) colloidal 
silica spheres of radius 285nm in deionized water, the 
following behavior was observed: For low nanoparticle 
concentrations the silica spheres aggregate, driven by 
the generic van der Waals attractions. At intermediate 
nanoparticle concentrations, the dispersion becomes sta- 
ble, whereas at higher concentrations the silica spheres 
aggregate again. The authors 3J attribute the initial sta- 
bilization to layering of the small nanoparticles near the 
colloidal surfaces. These "halos" occur because it is 
advantageous for the charged nanoparticles to be near 
the uncharged colloid surfaces. When two colloids then 
approach each other, their respective halos repel, pre- 
venting aggregation. The re-entrant aggregated phase, 
observed at higher nanoparticle concentrations, was at- 
tributed to normal entropic depletion attraction^. 

Clearly, a new route to stabilize colloids would have 
many potential applications. Indeed, nanoparticle halo- 
ing has already been used to enhance the self-assembly 



of 3-D colloidal crystals on patterned surfacesja]. Never- 
theless, even though this novel stabilization strategy has 
been demonstrated by experiment, many questions re- 
main about its generic applicability. To address these is- 
sues we carry out a systematic theoretical study of the ef- 
fective interaction (3 V b e /^ (r) 0| between colloids, induced 
by (much) smaller particles (P^—ksT is the reduced 
temperature.). We find fairly large regime of parameter 
space where (3V bb (r) is repulsive enough for stabiliza- 
tion, but this is usually followed by re-entrant attraction 
at higher small particle packing fraction. The picture 
that emerges is considerably more subtle than that of a 
simple static layer of adsorbed particles akin to steric sta- 
bilization. Instead, the nanoparticle halos are dilute, and 
in dynamic equilibrium with the bulk solution. Further- 
more, we observe no obvious change in their character 
when the re-entrant attraction kicks in, implying that 
this phenomenon is more complex than simple depletion 
attraction. 

The key quantity we study is the effective interaction^ 
@V b e b (r) between two spheres of diameter abb, induced 
by smaller spheres. Its properties are determined by the 
number density p s =N s /V of small particles and by the 
interactions /?$t, s (r) and /9$ ss (r); it is independent of 
the intrinsic interaction /?$;,{, (r) |]|. If $bb (f) is attractive 
and leads to aggregation, then introducing small particles 
that induce a (3V bb ^(r) repulsive enough to counteract 
(3&bb(r) will stabilize the colloids. 

The basic big-small and small-small interactions 
/3$jj(r) are modeled by a hard-core Yukawa (HCY) form 
which is versatile without having too many parameters 
to vary 7]. /3$y(r) = oo if r<a ij ; /3$y(r) = 0y(r) for 
r>Uij, where in each case r denotes the distance between 
the centers of the particles, and the Yukawa tail is 
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where o~bs = ^(o~bb+o~ss) with a ss the small particle hard- 
core diameter. By varying the size ratio q = a ss /ab s , 
packing fraction i] s = ^Trp s a^ s , and the 4 dimen- 
sionless potential parameters l3e ss , fttbs, A ss /cr ss and 
Abs/cfes, a wide variety of different physical situations 
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FIG. 1: (a): Comparison of simulations (from Q) and HNC 
calculations of the effective potential PV^ ! {r). The size- 
ratio q — 0.2, and the potential parameters X S a/<Jsa = 1/3, 
Abs/dss = 1/1.2, and /3e ss = 2.99 are kept constant, while 
/3e bs is varied, (b) HNC and DFtQ calculations for small 
size-ratios q. To facilitate comparisons, we plot Q/3V^/^ (r) 
and shift the curve for q — 0.1 up by +0.1. 

can be studied For repulsive interactions, for exam- 
ple, Eq. is a good model for charged suspensions 0- 

To restrict this vast parameter space somewhat, and 
inspired by the successful experiments Q, we choose 
e ss > and, initially, e bs = 0. The effective poten- 
tials @V bb (r) are calculated by using the two-component 
Ornstein Zernike (OZ) equations in the p b — > limit 
where they decouple, together with the hypernetted 
chain (HNC) integral equation closure^, leading to: 
pV b e b ff {r) = -p s J dr'h bs (r)c bs (\r' - r|), where h bs (r) = 
9bs( r ) — 1 with gbs( r ) the pair correlation function be- 
tween big and small particles, and c bs (r) is the direct cor- 
relation function, related to h bs (r) by the OZ equation[9j. 
In this limit, HNC has some important advantages |ll) 
over other popular integral equations such as Percus 
Yevick (PY) or Rogers Young (RY)0. For example, it 
is exact for the AO model|J| at all densities p s (PY is 
not Eft EH). Moreover, HNC is known to be particu- 
larly accurate for soft repulsive potentials of the type we 
are investigating 9| . To validate our method, we com- 
pare, in Fig. n the performance of HNC with several 
simulations [jj for e ss > 0, and find excellent agreement. 
Since we also want to study rather extreme size ratios, 
we compare, in Fig. to depletion potentials for hard 
spheres (HS) calculated with an accurate Density Func- 
tional Theory (DFT) approach [l^. Again HNC performs 
remarkably well. These results provide the confidence 
that, even if HNC is not perfectly quantitative, the trends 
we uncover will be robust, provided we limit ourselves 
to soft repulsions and low packing fractions 0|. Fortu- 
itously, this appears to be the regime where the nanopar- 
ticle haloing mechanism operates most effectively. 

In order to systematically investigate the conditions 
for which repulsive stabilization occurs, we calculated, 
with HNC, 0V bb (r) for a large number of parameter 



FIG. 2: Two typical examples of the effective potentials. 
First there is initial stabilization, and then a re-entrant at- 
traction. 

combinations. Two typical examples are shown in Fig. [2 
demonstrating the common pattern we find: for increas- 
ing packing fractions a maximum first appears close to 
contact, and continues to increase until at higher rj s a sec- 
ondary minimum appears that grows with rj s and rapidly 
moves to a separation of about one a ss . This sequence 
of initial stabilization followed by re-entrant attraction is 
similar to that seen in the experiments 0, and is found 
throughout the parameter regime we investigated. 

To further quantify the region of stability, we choose 
the following measure: For a given set /3e ss , X ss /a ss , and 
q, we calculate the effective potentials for different rj s , 
as done for Fig. The "stability window" is defined as 
A, ; = rf£ /rf s , where rj s s is the packing fraction above which 
the maximum of PV bb (r) is > 5 (leading to kinetic stabi- 
lization) , and 77" is the packing fraction below which the 
minimum of (3V bb ^(r) is < —2 (a conservative estimate 
of where short-range attractions induce aggregation 10]). 
The way our stability measure varies with potential 
parameters is depicted in Fig.0 from which some general 
trends can be extracted: The size of the window increases 
with increasing X ss /a ss and decreasing e ss and q. (In the 
HS limit (f3e ss — 0) we find no window of stability). For 
a number of points (A-D) in Fig. [3 we show the values of 
77" and 77". As expected, these packing fractions decrease 
with increasing (3e ss and \ ss / a ss since the small particles 
repel each other more and have a larger effective "size" . 

The effective repulsion is clearly related to the ac- 
cumulation of particles near the colloids (nano-particle 
halos0)- We define the "halo", as those particles be- 
tween r = a bs and r = r mmi the distance at which the 
pair correlation function g bs (r) has its first minimum. 
The number of particles Nhaio follows from integrating 
g bs {r) up to r min . The 2-D packing fraction is given by 
V 2D = ± ilT p 2D a 2 ss where p 2D = N halo /(A7ra 2 bs ). For all the 
parameters studied we find the same behavior depicted in 
the inset of Fig. OH The halo packing fraction rj 2D is linear 
with ?7 S and there is no change of slope or other obvious 
property marking either the beginning of stabilization at 



3 




FIG. 3: The ' 'ecjui-A^" lines denote different values of the 
stability window A,, = rff /77I for size-ratios q = 0.01 and q = 
0.05. For points A-D the stability window packing fractions, 
for q = 0.01, are listed in the format Inset: The 2D 

packing fraction rf D of a "halo" for different size-ratios q. The 
potential parameters are f3e as = 2, and \ ss /<j ss — 0.5. The 
symbols denote bulk packing fractions n" and rf£ for q = 0.01 
(circles), q — 0.05 (diamonds) and q = 0.1 (squares). 

r)g or re-entrant attraction at 77™. We have also inves- 
tigated other surface properties such as the adsorption 
r s = — J hbs(f)dr and the related surface tension t c [l5| . 
In contrast to polymeric depletants, where 7 S helps de- 
termine (r) . we observe no clear signatures of 
rjf or 77™ in the surface tension or the adsorption.. 

For small q we expect the gbs(r) to be similar within 
corrections 0(q^ 1 ), which explains why the rj 2D v.s. rj s 
curves are so close for different q. However, the stability 
windows, shown by the symbols in the inset of Fig. [31 dif- 
fer significantly: They are at lower rj s for smaller q, some- 
thing we observe more generally. This can be understood 
from an approximate Derjaguin0 argument valid in the 
small q limit. The potentials scale as 1/q times the force 
between two plates, and so stabilization (and re-entrant 
attraction) are achieved at lower packing fractions. How- 
ever, this doesn't easily explain why the window size also 
grows with decreasing q. 

The halos are very dilute at stabilization, and we have 
checked that all layer densities studied are well below that 
of any two-dimensional freezing transition. In fact, the 
layers shown in the inset of Fig. [3] are among the densest 
we investigated; for some 4> ss , rj 2D can easily be an order 
of magnitude lower at r?|. At these low packing fractions, 
the particles rapidly diffuse between halos and the bulk. 
In contrast to a steric stabilization mechanism, where 
the layers are static, we emphasize that this nanoparticle 
halo stabilization mechanism is dynamic. 

Further evidence against a naive picture of static layers 
comes from the re-entrant attraction. If the halos would 
become saturated, so that additional small particles can 
no longer segregate to the colloidal surface and instead 
act as depletants, then one might expect a linear depen- 
dence of the minimum of PV b /^(r) on r/ s as in AO0 



or HS01 depletion. Instead, the minimum in f3Vbb{r) 
grows initially as 77 J , closely resemb ling the behavior of 
the second minimum of HS systems |l2j, which suggests 
that both minima have a similar more complex origin in 
correlation effects. In fact, they would be directly re- 
lated if the potentials were interpreted in terms of a non- 
additive HS reference system with cr& s < \{obb + cr ss ), 
as explained in 0|. The repulsive effective interactions 
found in many other theoretical studies of f3V bb ^ (r) can 
also be qualitatively interpreted in this way (see e.g. Q 
and references therein for a discussion), suggesting that 
non-additivity ma y b e fruitfully used to interpret the re- 
entrant attraction 11]. 

One might argue that since adding an attractive 4>b s (r) 
should increase the number of particles in a layer, this 
should enhance the stabilization effect. However, we find 
more subtle scenarios. If we choose Afc s = A ss , to model 
residual charge on the large colloids, then for weak at- 
tractions the window indeed grows slightly. But, as Peb s 
becomes more negative, the potentials rapidly develop 
a large attractive component. This phenomenon, some- 
times called bridging for polymeric additivesQ], results 
from configurations where the two bigger colloids are 
both attracted to the same set of smaller particles 
An example of bridging is demonstrated in Fig.^Ja), and 
is representative of what we find more generally: the sta- 
bility window A, ; initially grows slightly, but then rapidly 
disappears, typically around (3eb s —1.5. 

On the other hand, a dramatic enhancement of the 
stabilization occurs for longer ranged colloid-nanoparticle 
attractions, as demonstrated in Fig. 0] for Xb s = 3A SS . 
The bridging effect is bypassed and the first minimum 
shifts up to positive absolute values, in fact for these 
parameters we find no re-entrant attraction within the 
range where we trust HNC0]. In general we find this 
effect for Xb s > X ss , but exactly where it kicks in depends 
q and the other potential parameters. In all cases studied, 
the 2D layer densities are still very low so that the "halos" 
are dilute; typically for more negative (3eb s bridging sets 
in againjllj. Of course when (3(f>b s (r) is attractive enough 
to induce static saturated layers, then the particles would 
be sterically stabilized. But technically this is a non- 
equilibrium effect: it only works if the colloidal particles 
are first isolated from each other on the timescale that 
the (saturated) layer forms. 

Finally, having systematically explored the effect of 
different parameters on f3V b /^ (r), we make some recom- 
mendations for experiments. For e/, s = the best sta- 
bilization should occur for modestly charged small ad- 
ditives (nanoparticles) since the stability windows are 
largest for small q, modest f3e ss and large \ ss /a ss . Small 
a ss are needed to enhance the latter parameter, since 
the (Debye) screening length is typically fixed by solu- 
tion conditions, e.g. X ss « 30nm for O.ItoM monovalent 
salt concentration. Another argument in favor of small 
particles concerns the dynamic nature of the layers. If the 
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FIG. 4: (a) Adding an attractive /3cj>b s (r) generates a deep 
minimum due to bridging effects when Xbs = A as . (b) A 
dramatic stabilization effect occurs for weak longer ranged 
attractions: At s = 3\ ss = Sa ss , and /3et s = —0.5. 

big particles are driven at velocities such that the smaller 
particles can no longer adiabatically follow, then the sta- 
bilization effect may disappear. Since the self-diffusion 
coefficients of the nanoparticles scale as 1/<t ss 0|, this 
again favors small particles. Moreover, larger windows 
A,; also enhance stability under halo fluctuations. 

Whereas adding charge to the colloids can destroy the 
stabilization effect, a modest but longer ranged attrac- 
tion P4>b s (r) can significantly enhance it. The latter ef- 
fect could be induced by residual van der Waals attrac- 
tions, although this recommendation must be tempered 
by the difficulty of adding van der Waals attractions 
between different species without simultaneously signifi- 
cantly increasing them between similar species. On the 
other hand, the advantage of small van der Waals attrac- 
tions (which are independent of Debye screening length) , 
is that the ratio Ah s /A ss can be tuned by changing salt 
concentration. This adds another handle for engineering 
effective potentials and concomitant phase behavior^]]. 

We observe the same general trends seen in the ex- 
periments of ref. Q , such as lower fy| with smaller q and 
values of t] 2D well below saturation. A direct quantitative 
comparison, however, is hampered by their use of polydis- 
perse small particles, and the difficulty in deriving accu- 
rate potentials (3(j>ij( r )- Taking their estimates (0e ss w 6, 
X SS /<7 SS « 0.6) we find a smaller window ps 2, at 
higher t] s than what they observed. The difference could 
stem from a small attractive (34>bs {f) or from polydisper- 
sity (preliminary calculations suggest that this lowers the 
effective rj^ and 77" (P. Bryk, private communication)). 

In conclusion, we discovered a substantial parameter 
regime where the addition of small (nano) particles can 
stabilize bigger colloids. Fortuitously, this occurs where 
the flexible HNC integral equation is most reliable. We 
usually find a stability window of packing fractions, above 
or below which the colloids aggregate again. The effects 



are significantly enhanced for weak longer ranged attrac- 
tive (34>bs{r). Although the stabilization is clearly related 
to the formation of diffusive accumulation layers around 
the bigger particles, we find no simple relationship to 
layer properties. This suggests these effects are related to 
more complex correlations. Colloidal stabilization by dy- 
namic nanoparticlc halos should be widely applicable and 
complimentary to existing steric and charge stabilization 
techniques This new mechanism may also have rel- 

evance to smaller scale biological interactions^^. [Note: 
upon completion of this work we became aware of a study 
by J. Liu and E. Luij ten, [cond-m at /0411278 which uses 
different techniques, but arrives at similar conclusions] 
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